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ABSTRACT 

We present a photometric study of the globular clusters in the giant elliptical galaxy M60 
in the Virgo cluster, based on deep, relatively wide field Washington CTi CCD images. The 
color-magnitude diagram reveals a significant population of globular clusters in M60, and a large 
number of young luminous clusters in NGC 4647, a small companion spiral galaxy north-west of 
M60. The color distribution of the globular clusters in M60 is clearly bimodal, with a blue peak at 
(C — Ti) = 1.37, and a red peak at (C — Ti) = 1.87. We derive two new transformation relations 
between the {C — Ti)q color and [Fe/H] using the data for the globular clusters in our Galaxy and 
M49. Using these relations we derive the metallicity distribution of the globular clusters in M60, 
which is also bimodal: a dominant metal-poor component with center at [Fe/H]= —1.2, and a 
weaker metal-rich component with center at [Fe/H]= —0.2. The radial number density profile 
of the globular clusters is more extended than that of the stellar halo, and the radial number 
density profile of the blue globular clusters is more extended than that of the red globular clusters. 
The number density maps of the globular clusters show that the spatial distribution of the blue 
globular clusters is roughly circular, while that of the red globular cluster is elongated similarly 
to that of the stellar halo. We estimate the total number of the globular clusters in M60 to be 
3600 ± 500, and the specific frequency to be Sn — 3.8 ± 0.4. The mean color of the bright blue 
globular clusters gets redder as they get brighter in both the inner and outer region of M60. This 
blue tilt is seen also in the outer region of M49, the brightest Virgo galaxy. Implications of these 
results are discussed. 
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1. Introduction 

Old globular clusters keep the fossil record for 
the early epoch of their host galaxies as well as 
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globular clusters themselves. By studying the 
age and metallicity of these globular clusters we 
can investigate the formation and early evolution 
of their host galaxies as well as globular clusters 
themselves. Globular clusters are distributed in a 
much wider region than the halo stars in their host 
galaxy, and thousands of them are found in giant 
elliptical galaxies (gEs). This feature, combined 
with easily-derived velocity, makes globular clus- 
ters powerful probes with which to study the struc- 
ture and kinematics in the outer halo of ne arby 
gEs as well as the i nner r egions fsee lLed ( 20031 ) and 
Brodie fc Strader ( 2006h . and references therein). 
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While the globular clusters in the inner re- 
gion of nearby gEs have been extensively stud- 
ied using the Hubble Spa ce Telescope (HST 

I2OOII 



leg using tne Mubbic bpa ce leiescope 
(iKundu fc Whitmorel bOOlt iLarsen et al. 



Peng et al. l2006l:lMies ke et al. IbOOdlStrader et al. I 
2006t [Harris et al. II2OO 6: Jord an et al. II2007I) . the 
globular clusters in the outer halo of nearby gEs 
have been studied using the wide f ield camera in 
the ground-based telesc op es (ILee fc Goislor 1993t 
McLaughlin et al. Ill99 3: Gcisle r. Lee, fc Kimll99a 



Lee. Kim fc Geisleil Il998t ,Ilhode fc Zepf 



2004: Dirsch et al. 120041: iHarris et al 



2001 



20041 : 



Tamura et al. I l2006al|bl : iBassino et al. I l2006h . 

However, the number of gEs for which the globular 
clusters were studied using the wide field camera 
is still small. 

Virgo, the nearest galaxy cluster, is one of the 
best targets for the study of globular cluster in 
gEs, because it includes several gEs that are al- 
most at the same distance from us. We have been 
carrying a long-term photometric study of glob- 
ular clusters i n Virgo gEs using the Washington 
filter system dLee fc Geisleil Il993t ICeisler I Il99l 
Lee. Kim fc Geisled Il998h . The Washinton sys- 
tem is known to be very sensitive to measuring 
the metallicity of the globular clusters and has a 
wide bandwidth so that it is ideal for studying 
the metallicity of the extra galactic globular clus- 



ters ( Geisler fc Forte 1990l ). The results for the 



globular clust ers in M87, th e centr al cD in Virgo 



were given in iLee fc Geisleil ([1993|) , and those for 



the globular clusters in M49, the brightest gEs in 
Virgo, were given i n jGeisler. Lee, fc KimI (ji996); 
Lee. Kim fc Geislei ( 19981 ). Here we present the 
results for the third Virgo gE, M60 in this series. 
M60 was selected because it is one of the brightest 
gEs in Virgo and is therefore expected to possess 
a rich globular cluster system. 

M60 (NGG 4649) is only shghtly less lumi- 
nous (A/y = —22.44 mag) than two brightest 
Virgo galaxies M87 (Afy = -22.62 mag) and M49 
(My = -22.83 mag). M60 is of morphological 
type E2, and is ultrav iolet-bright, emitting strong 
fiux at A < 25OOA (jBertola. Gapaccioli fc Oke 



19821 ). A recent Chadra image of M60 shows 
that the diffuse X-ray emission is detected out to 
about 3' fro m the center of M60 , with a circu- 
lar shape (,Humphrev et al. 20061 ). Basic infor- 



mation on M60 is listed in Table 1. We adopted 
a distance to M60, of 17.30 Mpc ((to - M)q = 



31.19 ±0.07) base d on the su r face b rightness fluc- 
tuation method in lMei et al. (2007) for which one 
arcsec corresponds to 84 pc. Foreground redden- 
i ng toward M60 is ver y small, E{B - V) = 0.026 
( Schlegel et al. I Il998l ). corresponding to E{C - 
Ti) = IM6E{B -V) ^ 0.051, yl(ri) = 0.071, 
and A{V) = 0.088. 

M60 has a companion SBc galaxy, NGG 4647, 
located 2. '5 from the center of M60 in the north- 
west direction (corresponding to a projected dis- 
tance of '^12.6 kpc for the adopted distance to 
M60). The radial velocity of NGG 4647 {v = 1422 
km s^^) is 305 km s~ ^ larger than that of M60 
(v = 1117 km s-^). IWhite. Keel, fc GoncehcT 



( 2000l ) found no foreground absorption due to 
NGG 4647 in the area of M60, and were unable 
to tell whether M60 or NGG 4647 is closer. 

Gouture et al. I ( 1991 ) performed the first pho- 



tometric study based on BV GGD imaging of the 
globular clusters in a small field (2'.1 x 3'. 4) of 
M60. They found a large dispersion in color dis- 
tribution, and a radial gradient in the mean cluster 
colors. They also found that the mean color of the 
globular clusters in M60 {{B - V) = 0.75) is 0.1 
mag redder than th at of M49 {{B ~ V) = 0.65). 



Harris et al. I (|l99lh determined the B band lumi- 
nosity function up to B ~26 mag of the globular 
clusters in this field, and found that the globular 
clusters follow a more extended spatial distribu- 
tion than the stellar light. Later, photometric 
studies based on HST/WFPG2 as well as ground- 
based images revealed that the M60 globular clus- 
ter system h as a clear bimodality in the color 
distri b ution (NeilsenI 19991 : Kundu fc Whitmore 



2001 



Larsen et al 
Forbes et al. I (|200 



2001; Forbes et al. 



l2004h . 



found from the analysis of 
wide field images obtained using Gemini/GMOS 
that the red globular clusters have a similar sur- 
face density distribution to that of the stellar light, 
which is steeper than that of the blue globular 
clusters. In addition, they derived a value for the 
specific frequency of the M60 globular clusters of 
Sn ~ 4.1 ± 1.0, which is lower than the value 
Sn ~ 6.7 given in lAshman fc Zepf I (jl998l ) based 
on inferior data. 

Recently, it has been found that M60 is one 
of the galaxies that show a "blue-tilt" in the 
color-magnitude diagram of its globular clus- 
ters, in which the brighter the blue globu- 
lar clusters are, the redder they are in the 
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mean (iHarris et al I l200d IStrader et al. I l2006t 



Mieske e t al. I I2006D. iSarazin et al. I (|2003h and 
Randall. Sarazin fc IrwinT ^ 20041 ) detected some 
discrete source s in M60 in the X-ray ban d usin, 



Chandra, and iRandall. Sarazin fc Irwin" ( 200 



found that roughly 47% of the X-ray discrete 
sources is identified with globular clusters. By 
cross-correlating Chandra poin t sources and opti- 



cal globular cluster candidates. iKim et al. ( 20061 ) 



found that the mean probability for a globular 
cluster to harbor a low mass X-ray binary (LMXB) 
in M60 is about 6.1±1.0 %, and that this probabil- 
ity for the red GCs is much larger th an that for the 
blue g lobular clusters. In addition, IPierce et al. 
pood ) pubhshed a spectroscopic study of 38 glob- 



ular clusters in M60 b ased on data ob t ained using 
Gemini/GMOS, and [Bridges et afl (l2006l) pre- 
sented a study of the globular cluster kinematics 
of M60 using these data. 

In this paper we present a photometric study 
of globular clusters in M60 using deep, relatively 
wide field CCD photometry. This study covers a 
field of M60 that is much larger than any previous 
studies on globular clusters in M60, and is using 
the Washington system that has a better sensitiv- 
ity for measuring the metallicity compared with 
most other systems. This study supplements also 
our k inematic studies of the globular clus ters in 



M60 (|Lee et al. 112001 iHwang et al. Il2008h . This 



paper is organized as follows. In Section 2 we de- 
scribe our observations and data reduction. In §3 
we present the color-magnitude diagram and color 
distribution of the globular clusters, and compare 
the surface photometry of the stellar halo with the 
structure of the globular cluster system. We also 
investigate the radial variation of the mean mag- 
nitudes and colors of the globular clusters, and es- 
timate the total number and specific frequency of 
the globular clusters. In §4 we discuss our results 
and their implication in comparison with other 
studies. Primary results are summarized in the 
final section. 

2. Observations and Data Reduction 

2.1. Observations 

CCD images of M60 were obtained on the pho- 
tometric nights of April 9 and 10, 1997 (UT) at 
the KPNO 4 m telescope using the 2048 x 2048 
pixels CCD camera at the prime focus. Washing- 



ton C and Kron-Cousins R filters were used. The 
Kron-Cousins R filter has a very similar effective 
wavelength to that of the Washington Ti filter, 
but with a much wider bandwidth re s ulting in 
three times the sensitivity ( Geisler |[l996[) . Geisler 



(Il996l) showed that Ti reproduces R magnitude 
very well: R = Ti + 0.003 - 0.017(C - Ti) with 
rms=0.02. So we used the R filter as an alterna- 
tive to the Ti fi lter, as done in most recent Wash- 
ingto n studies dCeisler 19961 : Lee. Kim fc Geisler 



1998 : 



Bassino et al 



Dirsch et al 
2n06h 



20041: 



Harris et al. 



20041: 



The observation log is given in Table 2. The size 
of the field of view is 16'. 4 x 16'. 4, and the pixel 
scale is 0.47 arcsec pixel"^. Exposure times are 
100 s and 4 X 1500 s for C, and 60 s and 3 x 1000 s 
for Ti. Each long exposure image was taken with 
dithering of 10 to 28 arcsec. The seeing ranged 
from 1.1 to 1.5 arcsec. I n additio n , seve ral Wash- 
ington standard fields in I Geisler" ( 1996[ ) were ob- 
served during the observing run. 

2.2. Point Source Photometry 

Each frame was trimmed, bias-subtracted and 
flat-fielded with twilight skyflats for C and dome 
flats for Ti using the IRAF software. The indi- 
vidual long exposure images in each filter were 
then shifted to a common center and medianed 
together. Figure 1 displays a grayscale map of the 
short exposure Ti image of M60. 

Globular clusters at the distance of M60 ap- 
pear as point sources in our images. First we 
subtracted the stellar halo of M60 from each of 
the original images for better detection of the 
sources as follows. We created a model image of 
M60 using the ELLIPSE task in IRAF/STSDAS. 
Then we subtracted the model image from the 
original image to remove the stellar halo light, 
and then added a constant equal to the mean 
value of the sky background in the original im- 
age. We used the resulting images for the de- 
tection of objects with the digital photome- 
try p rogram DAOPHOT/ALLFRAME (.Stetsonl 
19941 ). We used threshold=3cr, sharpness=0.2- 
1.0, roundness=-0.7 to 0.7 for detection in the 
KPNO images. Instrumental magnitudes of the 
detected objects in the images were derived using 
the point-spread function (PSF)-fitting routine 
in DAOPHOT/ALLFRAME. To take care of the 
variation of the PSF we adopted a quadratically 
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variable PSF, and applied the aperture correc- 
tion to the PSF-fitting magnitudes. The value 
of the aperture correction was derived from the 
difference between the aperture magnitudes and 
PSF-fitting magnitudes for isolated bright stars: 
0.014 and 0.008 for short and long exposure C im- 
ages, respectively, and 0.033 and -0.009 for short 
and long exposure Ti images, respectively. 

The Ust of detected objects includes both point 
sources and extended sources. We selected the 
point sources among the detected sources using 
the morphological classifier r„2 moment defined as 

r„2 = - sky)/^,[{h - sky)/irf + 0.5)]f\ 

where fj is the radial distance of the i-th pixel 
from the center of the source, li is the intensity 
value at the i-th pixel, and sky is the local back- 
ground value derived from a more distant annulus 
(|Kron|[l98Q) . The FWHMs of the point sources 
vary approximately as a quadratic function of the 
radius from the center of the image. So r_2 was 
normalized to the value for the center of the im- 
age. From the artificial star experiment to be de- 
scribed below, we decided to consider the sources 
with r_2 < 1.27 to be point sources. Misclassi- 
fication rates are negligible for the bright objects 
with Ti < 23 and increase with magnitude, reach- 
ing about 20% at Ti = 24. The central regions at 
r < 1 arcmin were saturated in the long exposure 
images so that we used the photometry of the ob- 
jects in the central region derived from the short 
exposure images. 

For this reason, the KPNO images are of lim- 
ited use for the study of the globular clusters in 
the central region of M60. Therefore we analyzed 
F555W and F814W images of the central region 
(2' X 2') of M60 in the HST/WFPC2 archive to 
supplement the KPNO data. Hereafter we call V 
and / for F555W and F814W, respectively. The 
position of t he HST/WFPC2 field is marked in 



Figure 1 (see Kim et al. ( 2006h for the journal of 



its observation) . We created a model image of M60 
using the ELLIPSE task in IRAF/STSDAS. Then 
we subtracted the model image from the original 
image to remove the stellar halo light. We de- 
tected objects in the WFPC2 images and derived 
the photometry of the detected objects using the 
digita l photometry software HSTPHOT (|Dolphinl 
2000l ). We used 3.5 a as a threshold for detec- 
tion in the HST images. We used the radii of the 
aperture of 3 pixels for PC chips and 2 pixels for 



WF chips to get the aperture magnitudes of the 
detected objec ts, and applied the aperture correc- 
tion given by iKundu &: Whitmord ( 200ll ) to get 
the total magnitudes of the detected objects. We 
selected the starlike sources in the list of the ob- 
jects returned by HSTPHOT using r_2. We used 
the HST photometry for the analysis of the central 
region (r < 1'.5) on one side of M60 for which the 
KPNO photometry is poor. We applied the same 
procedure to the images of another HST/WFPC2 
field, at 5' north from the center of M60 (called the 
nort h HS T/WFPC2 field), as marked in Figure 1 
(seeE im et al. I ()2006l ) for the basic information of 
this field). It is found that the number of globular 
clusters in this field is so small that they are of 
limited use for this study. 

2.3. Standard Calibration 

We derived the transformation equations for 
standard calibration from the photometry of the 
Washington standard stars observed during the 
same night. We obtained the instrumental mag- 
nitudes of the standard star s using the ap erture 



radius of 7.5 arcsec as used in lGeislerl(|l996h . The 



standard transformation equations we derived are: 
(April 9) Ti=ti + 0.034(c - h) - 0.122A: -I- 0.282 
with rms=0.021 and N=42, and (C - Ti) = 
1.062(c - ti) - 0.263X - 0.405 with rms=0.018 
and N=44, and (April 10 ) Ti = ii -I- 0.034(c - 
h) - 0.122X + 0.269 with rms=0.028 and N=39, 
and (C-Ti) = 1.062(c-ti)-0.263X-0.425 with 
rms=0.030 and N=39, where the upper case let- 
ters represent the standard magnitudes, the low- 
ercase letters the instrumental magnitudes (with a 
DAOPHOT system zero point of 25.0), and X the 
air mass. We transformed the instrumental mag- 
nitudes of the sources onto the standard system 
using these transformation equations. We selected 
and listed in Table 3 the CTi photometry of 4497 
point sources with a{C — Ti) < 0.3 measured in 
the KPNO images, and Figure 2 displays the mean 
errors of Ti and (C — Ti) versus Ti magnitude. 

2.4. Completeness of the Photometry 

We estimated the completeness of the KPNO 
photometry using DAOPHOT/ ADDSTAR that 
was designed for the artificial star experiment. We 
generated a set of artificial stars for which the 
color-magnitude diagrams and luminosity func- 
tions are similar to the observational one, using 



4 



the PSFs derived from the long exposure real im- 
ages. We did not include the central region at 
r < 1' that was saturated in the real images. Then 
we added them to the real image avoiding the po- 
sition of detected real objects. We added 1200 
artificial stars to each pair of C and Ti images in 
a set of 50 pairs so that the total number of added 
artificial stars is 60,000. Then we applied the same 
procedure of photometry to the artificial images as 
used for the real images, and estimated as the com- 
pleteness factor, the number ratio of the recovered 
artificial stars and the added artificial stars. 

Figure 3 displays the completeness for C and Ti 
we derived for the point sources. We plotted the 
completeness for the entire region and five bins 
in radial distance: 1' < r < 2', 2' < r < 3', 
3' < r < 4', 4' < r < 5', and 5' < r < 6'. It is 
seen that the completeness is higher than 90 % for 
Ti = 23 (C = 24.3) for the outer region at r > 2' 
and Ti = 22.8 (C 24.0) for the inner region at 
V < r < 2'. The photometric limit levels for 50 
% completeness are Ti = 24.4 (C = 25.4) for the 
outer region at r > 2' and Ti = 23.9 (C = 25.0) 
for the inner region at 1' < r < 2'. The com- 
pleteness varies little depending on radius for the 
outer region at r > 2', while it is somewhat lower 
for the inner region at r < 2' compared with the 
outer region. We derived the mean photometric 
errors from the difference between the magnitudes 
of the added objects and recovered objects, finding 
that they are similar to the measured photomet- 
ric errors for the real objects as seen in Figure 2. 
We derived similarly the completeness for V and / 
using HSTPHOT, which is plotted also in Figure 
3. The completeness for the HST photometry is 
higher than to 95 % for V < 23.5 mag for all the 
radial range of the HST data. 

2.5. Comparison with Previous Photome- 
try 

There is no previous Washington photometry 
of the sources in M60. However there are a few 
photometry of the sources in M60 using the dif- 
ferent filter systems in the literature. We derive 
the transformation relations between our Wash- 
in gton photometry and the gi photometry given 



bv lForbes et al. I (|2004[ ). 



iForbes et al. 



(j20041 ) presented gi photometry 
of the globular clusters in a field of about 90 
square arcmin including M60, based on the CCD 



images taken using the Gemini North telescope. 
Our field of view covering 256 squ are arcmin is 
about 2.8 times larger than that of Forbes et al. I 
(|2004l ). We compared (C - Ti) colors and {g - i) 
colors for 396 objects with Ti < 22 mag in com - 
mon between this study and lForbes et al.1 ( 20041 ). 
as displayed in Figure 4. Figure 4 shows that the 
relation between {g — i) colors and (C — Ti) col- 
ors can be fit by a triple linear relation as follows: 
{g-i) = 0.735(C-Ti) -0.121 with rms^O.llO for 
0.2 < (C-Ti) < 1.2, [g-i) = 0.527(C-ri)+0.119 
with rms=0.092 for 1.2 < (C - Ti) < 2.4, and 
(5 - i) = 1.372(C - Ti) - 1.943 with rms=0.243 
for 2.4 < (C - Ti) < 4. In addition the trans- 
formations between magnitudes are derived to be: 
g = C+0. 340-0. 548(C-ri) with rms=0.067, and 
i^Ti+ 0.429 - 0.199(C - Ti) with rms=0.070. 

We also derived the transformation between 
(C — Ti) colors and {V — I) colors for the objects 
with Ti < 22 mag in common between the KPNO 
images and HST/WFPC2 images, as shown in Fig- 
ure 4. They show a linear relation for the color 
range of 0.9 < (C-Ti) < 2.4. Linear fitting to the 
data for 72 objects with small errors (cr(C — Ti) < 
0.1 ) yields {V - I) = 0.407(C - Ti) + 0.459 with 
rms=0.047 for 0.9 < (C-Ti) < 2.4. This is in very 
good agreement with the result deri ved from the 
photo metry of M49 (NGC 4472) by iLee fc Kim 
(|200n[ ). (V ~ I) = 0.443(C - Ti) + 0.396 for 
(C — Ti) < 2.6. We also derived an approximate 
relation for the magnitudes, V — Ti — 0.106 + 
0.243(C-Ti) with rms=0.100 for 66 objects with 
0.9 < (C - Ti) < 2.4. 

3. Results 

3.1. Color-Magnitude Diagram 

Figure 5 displays the color-magnitude diagrams 
of the measured point sources in the KPNO im- 
ages of M60 (in three regions at 1'.5 < r < 7', 
7' < r < 9', and 9' < r < 10' where r is the pro- 
jected galactocentric distance) as well as those in 
the HST/WFPC2 images (for r < 1'.5). We ex- 
cluded the objects within a circular region of radial 
distance of 1' from the center of NGC 4647 from 
the KPNO data in Figure 5. We transformed the 
{V - I) colors of the objects in the HST WFPC2 
images into (C — Ti) colors using the transforma- 
tion equation in the previous section. 

Three general kinds of objects are seen in Fig- 
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ure 5: (a) globular clusters of M60 residing in 
the broad vertical feature in the color range of 
1.0 < (C - Ti) < 2.4, (b) a small number of 
bright foreground stars at (C — Ti) < 0.9 and 
(C — Ti) > 2.6, and (c) faint blue unresolved back- 
ground galaxies with Ti > 23 mag. Two vertical 
features are visible in Figure 5(b) at 1'.5 < r < 7': 
blue globular clusters (BGCs) and red globular 
clusters (RGCs). We selected a sample of bright 
globular clusters with 1.0 < (C - Ti) < 2.4 and 
19 < Ti < 23 (-12.26 < Mr, < -8.26) for the 
analysis, considering the following points: 1) the 
color range of the known globular clusters in other 
galaxies is (C — Ti)o « 1 to 2.4, and the fore- 
ground reddening for M60 is small; 2) the peak 
luminosity of the globular clusters is estimated to 
be somewhat fainter than T sa 23 at the distance 
of M60 (see also Section 3.7); 3) the mean photo- 
metric error is smaller than cr(Ti) — 0.05 and the 
incompleteness in our photometry is estimated to 
be minor for Ti < 23; and 4) the contamination 
due to background galaxies is estimated to be very 
small for Ti < 23. The sample of globular clus- 
ters is separated around (C — Ti) = 1.7, since the 
color boundary between the two components is es- 
timated to be at (C — Ti) = 1.7 as described in the 
following section. So we divided the entire sam- 
ple of globular clusters into two classes: the BGCs 
with 1.0 < (C - Ti) < 1.7 and the RGCs with 
1.7 < (C - Ti) < 2.4. A small number of the 
BGCs are still seen even at 9' < r < 10', while 
very few RGCs are visible in the same region. 

Figure 6 displays the color-magnitude diagrams 
of the point sources within a radius of one arcmin 
from the center of the companion spiral galaxy 
NGC 4647 (upper panel) and in a "control field" 
with the same area (lower panel). The region for 
NGC 4647 is about 2.5' from the center of M60 so 
that it is unaffected by the saturation of the M60 
nucleus in the image. The control field is located 
at the same radial distance from the center of M60, 
but in the opposite direction. We also checked 
the color-magnitude diagrams for two other con- 
trol fields at the same radial distance but at dif- 
ferent position angle (45 degrees and 315 degrees, 
respectively). They are very similar to the color- 
magnitude diagram for the chosen control field in 
Figure 6, so that the chosen control field can be 
considered to represent a control field for the sta- 
tistical analysis. The number of objects within 



the boundary of bright globular clusters in the 
field of NGC 4647 (marked by the two rectangles) 
is slightly smaller than that in the control field, 
showing that most of the objects within the rect- 
angles are probably globular clusters belonging to 
M60. Most of the point sources in NGC 4647 are 
much bluer than globular clusters, but as bright 
as globular clusters of M60. 

There are seen many more sources that are 
slightly extended, compared with the point sources, 
in NGC 4647. We selected the extended sources 
using the criterion 1.27 < r_2 < 2 (note r_2 < 
1.27 for the point sources). They are slightly 
larger than the point sources and were included in 
the output fist of the DAOPHOT/ALLFRAME. 
They have larger values for x and sharpness than 
the point sources, indicating that their PSF-fitting 
magnitudes are less reliable than those of the point 
sources. Although their r_2 values show that they 
are slightly larger than the point sources, it is not 
easy to estimate their size in our images. We de- 
rived the photometry of these objects, plotting 
them in Figure 6. Most of these slightly extended 
sources are also much bluer than globular clus- 
ters, but as bright as globular clusters of M60. 
These blue objects are not seen in the control 
field. Therefore these point sources and slightly 
extended sources with blue colors are considered 
to be young luminous star clusters or associations 
in NGC 4647. 

Recent studies show several evidences of inter- 
action between N G C 464 7 and M60, although 
Sandage fc Bedke I (|l994l ) described that the 
two galaxies are not interacting: (1) The mor- 
phology of NGC 464 7 is clearly asymmetric 



(jKoopman et al. I 120011 ): (2) The inner region 



of M60 has a strong rotational support in con- 
trast to other giant elliptica l galaxies, and has an 



asymmetric rotation c urve (jPinknev et al. 1 120031: 



de Bruvne et al. Il200lh : and (3) there is seen a fil- 
ament extending to the northeastern edge of M60 



in th e X-ray image (jRandall. Sarazin fc Irwin 
20061 ). Therefore we suggest that the young clus- 



ters in NGC4647 might have formed during a re- 
cent interaction of NGC 4647 with M60. 

3.2. Color Distribution of the Globular 
Clusters 

Figure 7 displays the color distribution of the 
bright point sources with 19 < Ti < 23 at 1' < 
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r <T (most of which are globular clusters) in the 
KPNO images. We did not include the objects 
at r < 1' because they suffer from incompleteness 
problem. We also derived the color distribution of 
the background objects with the same magnitude 
range at 9' < r < 10', finding the contribution of 
these background objects is negligible as shown in 
Figure 7. 

The color distribution is clearly bimodal. KMM 
tests of the data show that the probability that 
the color distribution is bimodal over unimodal is 
higher than 99.9 %. We determined the follow- 
ing parameters for the best-fit double Gaussian 
curves describing the color distribution, using a 
maximum-likelihood met hod through the KMM 
mixtu re modeling routine ( Ashman. Bird, fc Zepf I 
1994 ): a primary component with center at (C — 
Ti) = 1.37 and width a — 0.16 and a secondary 
component with center at (C — Ti) = 1.87 and 
a = 0.23. We did not assume equal dispersions 
for the blue and red globular subpopulations in 
these KMM fits and they do indeed appear to be 
significantly different. 

The minimum between the two components is 
found to be at (C — Ti) = 1.7, which was used 
for dividing the entire sample into BGCs and 
RGCs. This boundar y color is slightly redder 
than that used for M49 iGeisler. Lee, fc Kimlll996l : 



Lee. Kim fc Geisle r 1998), but is consistent with 
Couture et al. I ([1991,) 's finding based on BV pho- 
tometry of globular clusters in M60. They also 
found that the mean color of the globular clusters 
in M60 {{B -V) = 0.75) is 0.1 redder than that 
of M49 {{B - y) = 0.65). 

We investigate the radial variation of the color 
distribution of the bright globular clusters after 
subtracting the background level derived from the 
region at 9' < r < 10' in Figure 8. In the 
central region at r < 75", the RGC component 
looks stronger than the BGC component, while 
the BGC component begins to dominate for r > 2' 
and steadily strengthens its domination in the 
outer region. The RGC component is barely seen 
in the bacgkround region at 9' < r < 10'. 

We also derived the variation of the ratio of 
the number of the bright BGCs to the number of 
bright RGCs with r < 7', plotting it in Figure 
9. Figure 9 shows that this ratio keeps increas- 
ing as the galactocentric radii increase, and is fit 
well by the hnear relation N [BGC) /N {RGC) = 



0.43(±0.04)r -I- 0.47. This ratio becomes unity at 
r = 74". 

3.3. Galaxy Surface Photometry 

We derived the surface photometry of M60 from 
the KPNO images using the ELLIPSE task in 
IRAF/STSDAS. First we masked out bright fore- 
ground stars and nearby galaxies including NGC 
4647 in the original images. Then we applied EL- 
LIPSE for ellipse fitting of the isophotes of M60 
to the resulting images. We used the very outer 
region at r ~ 10' in the corner of the image to 
estimate the background level. We measured the 
median intensity value for each of three regions in 
the corner and took the median of the three mea- 
dian values as the background level. The central 
part (r < 3 arcsec) of the short exposure Ti image 
of M60 was saturated so that we could not derive 
the color for this region. The radial profiles of the 
surface color were obtained using the same struc- 
tural parameters for both C and Ti, which were 
derived from the Ti images. The errors for the 
surface brightness magnitudes are those that are 
given by ELLIPSE, not including the error for the 
background estimate. 

Figure 10 displays the radial profiles of the 
surface brightness magnitudes, (C — Ti) color, 
ellipticity and position angle (PA) of M60 as a 
function of major radius rmaj- We fit the surface 
brightness profiles with a de Vaucouleurs r^/^ law 



for the range of 3". 8 < 



maj 



< 410".l, using linear 



least squares fitting: fi{C) = 2.568(±0.011)ri/4 + 
14.955(±0.032) with rms=0.065, and ^j.{Ti) = 
2.649(±0.012)ri/4-hl2.832(±0.034) with rms=0.069, 
which are also plotted by the solid lines in the 
same figure. It is found that the surface bright- 
ness profiles are well fit by a de Vaucouleurs r^^^ 
law. The effective radius and standard radius for 
C-band (where ^(C) = 25.0 mag arcsec"^) are de- 
rived to be 97 arcsec and 242 arcsec, respectively, 
corresponding to linear sizes of 8.15 kpc and 20.3 
kpc. 

The (C — Ti) color gets bluer in the inner region 
(r < 150") as the galactocentric radius increases, 
and remains almost constant in the outer region. 
This trend for the inner region is consistent with 
the radia l variation of (U — R) and {B — R) colors 
given by IPeletier et al. ( 1990) as seen in Figure 
10. 
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The color trend at 3". 2 < r < 151". 4 can 
be fit by the log- linear relation (C — Ti) = 
-0.162(±0.011)logr+2.160(±0.032) with rms=0.065. 
The ellipticity increases rapidly from 0.1 to 0.2 in 
the central region, and stays at an almost con- 
tant value in the outer region. The values of the 
ellipticity at the effective radius and standard ra- 
dius are derived to be, respectively, 0.21 and 0.22. 
The PA increases from 90 deg to 100 deg in the 
central region, and changes slowly in the outer re- 
gion. The values of PA at the effective radius and 
standard radius are, 106 and 110 deg, respectively. 

3.4. Spatial Distribution of the Globular 
Clusters 



We have investigated the spatial structure of 
the globular cluster system in M60. Figure 11 dis- 
plays the spatial distribution of the bright globular 
clusters with 19 < Ti < 23 mag (all globular clus- 
ters, BGCs and RGCs) as well as foreground stars 
found in the KPNO images. We also plotted the 
very blue bright clusters (YC) with (C-Ti) < 0.5 
and 19 < Ti < 23 mag. We selected the red 
point sources with (C - Ti) > 2.5 and Ti < 23 
mag as foreground stars for comparison. The spa- 
tial distribution of these objects is seen to be uni- 
form in Figure 11(b), showing indeed that they 
do not belong to M60. We also created number 
density maps of the globular clusters, which were 
smoothed and displayed in Figure 12. In Figure 
12(b) we also display a grayscale map of the short 
exposure Ti image of M60 that shows the distribu- 
tion of stellar light, which is to be compared with 
the globular clusters. 

Several notable features are seen in Figures 11 
and 12. First, the spatial distribution of all glob- 
ular clusters is roughly circular, and it shows a 
strong central concentration. Second, the spatial 
distribution of the BGCs is roughly circular and it 
extends farther than that of the RGCs. Third, the 
spatial distribution of the RGCs is somewhat elon- 
gated along the east-west direction, which is con- 
sistent with the position angle of the halo of M60. 
We derive e(ellipticity)= 0.05 ± 0.02, 0.03 ± 0.02, 
and 0.09 ± 0.04, respectively, for all globular clus- 
ters, BGCs and RGCs at 1'. 5 < r < 7^ using the 
disper sion ellipse method of iTrumpler fc Weaver 
(jl953l ). Fourth, very blue, presumably young, 
bright clusters are located mostly within one ar- 
cmin from the center of NGC 4647 (inside the 



small circle at the position of NGC 4647 in Figure 
11). 

The Chanda X-r a y ima ge of M60 (in Fig. 1 
of Humphrev et al. ( 20061 )) shows that the X-ray 
emission of M60 is almost circular. Therefore the 
global structure of M60 seen in the X-ray image 
is closer to the spatial distribution of the BGCs 
rather than to that of the RGCs. This indicates 
that both the BGCs and the X-ray emitting hot 
gas may trace the dark mass of their host galaxies. 

3.5. Surface Number Density Profiles of 
the Globular Clusters 

We have derived the radial profiles of the sur- 
face number density of bright globular clusters 
with 19 < Ti < 23 mag, using the KPNO data 
for the outer region at r > 1'.5 and the HST data 
for the central region at r < 1'.5. First we checked 
the radial profiles of mean counts per unit area for 
the objects with the same range of magnitude and 
color as the globular clusters in the KPNO im- 
ages, finding that they get almost flat at r > 8'. 
Therefore we derived the background levels from 
the mean surface number density of these objects 
at 9' — 10': 1.988 ± 0.363 per square arcmin for all 
globular clusters, 1.790 ± 0.344 per square arcmin 
for the BGCs, and 0.199±0.115 per square arcmin 
for the RGCs. Then we subtracted these back- 
ground values from the original number counts to 
produce the radial profiles of the net surface num- 
ber density of globular clusters, which are listed 
in Table 4. We derived the surface number den- 
sity profiles in the core region using the HST data 
with the same criteria as those for the KPNO data. 
We estimated the background contribution for the 
HST data using the result derived from the KPNO 
data, finding that it is negligible. 

Figure 13 displays the radial profiles of the sur- 
face number density of globular clusters in com- 
parison with the C-band surface brightness profile 
of M60. Several features are of note in Figure 
13. First, the surface number density profile of 
all globular clusters extends farther than the sur- 
face brightness profile of the galaxy halo. Second, 
the surface number density profile of the RGCs 
agrees approximately with the surface brightness 
profile of the galaxy halo, although the halo is 
steeper. Third, the surface number density pro- 
file of the BGCs extends farther than that of the 
RGCs. Fourth, the surface number density profiles 
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of the globular clusters in the outer parts at 0'.5 < 
r < 1' are fit approximately both by the de Vau- 
couleurs r^/-* law: logti = -1.974(±0.047)ri/4 + 
3.571(±0.051) for all globular clusters, logcr = 
-1.717(±0.068)ri/4+2.984(±0.074) for the BGCs, 
and logcr = -2.262(±0.050)ri/4 + 3.580(±0.052) 
for the RGCs; and by a power law: logcr = 
-1.285(±0.034) logr + 1.563(±0.011) for all glob- 
ular clusters, logcr = -1.128(±0.040) logr + 
1.237(±0.014)forthe BGCs, and logcr = -1.452(±0. 
1.281(±0.016) for the RGCs. Effective radii are 
8'.09, 14'. 14, and 4'. 70 for all globular clusters, 
BGCs and RGCs, respectively. Thus the sur- 
face number density profile of the RGCs is much 
steeper than that of the BGCs. 

Finally, the surface number density profiles of 
the globular clusters are flat in the central re- 
gion at r < 0'.5. This result is affected little by 
the incompleteness of our photometry, because the 
completeness is higher than 95 % for F < 23.5 
mag for all the radial range of the HST data. 
This flattening indicates the globular cluster sys- 
tem is dynamically relaxed in the central region 
in the sense that the radial profile of th e cen- 



tral r egion can be fit by the King model (jKing 
19621 ). Figure 14 shows the results derived from 
fitting the data for r < 2' with a King model. 
From error-weighted fitting we derive core radius 
Tc = 0.84' and concentration parameter c = 3.91 
for all globular clusters, Tc = 1.06' and c = 3.92 
for the BGCs, and = 0.71' and c = 3.57 for 
the RGCs. From equal-weighted fitting we derive 
similar results: Tc = 0.73' and c = 2.83 for all 
globular clusters, rc — 0.86' and c — 2.65 for the 
BGCs, and = 0.66' and c = 2.98 for the RGCs. 
However, it is not easy to determine reliably the 
values of both parameters. If we adopt a fixed 
v alue, c = 2. 5 , as d one for the case of M87 by 



Kundu et al. 



(1993), we derive = 0.85' for all 
globular clusters, Tc = 1.06' for the BGCs, and 
Tc = 0.72' for the RGCs. Thus the core radius 
for the RGCs is much smaller than that for the 
BGCs. It is also noted that the surface number 
density profile of the RGCs at the outer area is 
reasonably fit by the King model, while that of 
the BGCs at r > 3' shows a clear excess over the 
King model. The c entral flatten ing is also seen 
in the case of M49 (|Lee fc Kim|[2000j and M87 
(|Kundu et al. II1999I) both of which were also based 
on the HST/WFPC2 data. The flattening for the 



globular clusters can be explained by the intrinsic 
property of the glo bular cluster formation epoch 
( Harris et al. 1998f ) or by the effect of the globu- 
lar clus ter destruction du ri ng the dynamica.1 evo - 
lution (jCote et al. Ill998l ). iKundu et ahl (|l999[ ) 
concluded there is little evidence supporting the 
latter in the case of M87, noting that the luminos- 
ity function of the globular clusters in M87 does 
not vary spatially within about 1'.5 from the cen- 
050) Idsgr-^^f M87. Therefore the flattening for the M60 
globular clusters may be the relic of the formation 
epoch of the globular clusters. 

3.6. Radial Variation of Mean Magnitudes 
and Colors 

Figure 15 displays the colors of the globu- 
lar clusters with 19 < Ti < 23 versus galac- 
tocentric distance derived from the KPNO data 
only. Mean and median colors of the globu- 
lar clusters in each radial bin are represented by 
the larger symbols. The mean color of the en- 
tire sample of globular clusters shows a clear ra- 
dial gradient, while those of the BGC and RGC 
show little, if any, radial gradients. Linear least 
squares fitting for the range of 1' < r < 9' 
yields (C - Ti) = -0.033(±0.003)r + 1.718 with 
rms=0.020 for aU globular clusters, (C - Ti) = 
-0.015(±0.003)r-h 1.457 with rms=0.012 for the 
BGCs, and (C - Ti) = -0.008(±0.002)r -h 1.984 
with rms=0.009 for the RGCs. The samples at 
7- < 1' in the KPNO data suffer from incomplete- 
ness so that they were not included for fitting. The 
mean color of the stellar halo is much closer to that 
of the RGCs than to that of the BGCs for r < 8', 
but is 0.2 to 0.3 mag bluer than that of the RGC. 
This is in contrast to the case of M49 where the 
color profiles of t he RGCs agree well with that of 
the steUar halo (|Lee. Kim fc Geisleij [l998l ) . The 
mean color for the RGCs depends on the color 
range adopted for taking a mean so that the color 
difference between the halo and the RGCs is not 
considered to be significant. 

Figure 16 displays the C and Ti magnitudes of 
the globular clusters with 19 < Ti < 23 mag ver- 
sus galactocentric distance from the KPNO data 
only. Mean magnitudes of the globular clusters in 
each radial bin are represented by the larger sym- 
bols. The mean magnitudes of the globular clus- 
ters with r > 1' vary little as a function of galacto- 
centric distance. The upper (bright-side) envelope 
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of the globular cluster magnitude distribution also 
does not show any clear systematic radial gradi- 
ent. Linear least squares fitting for the range of 
1' < r < 9' yields Ti = 0.007(±0.008)r + 21.945 
with rms=0.053 for all globular clusters, Ti = 
0.012(±0.017)r + 21.881 with rms=0.058 for the 
BGCs, and Ti = 0.033(±0.013)r + 21.866 with 
rms=0.050 for the RGCs. The samples at r < 1' 
suffer from incompleteness problem so that they 
were not included for fitting. Thus we find little 
or no dependence of the mean magnitude of the 
globular clusters on galactocentric distance. 

3.7. Luminosity Function and Specific Fre- 
quency 

Figure 17 displays the luminosity function of 
the globular clusters in M60. We plot the luminos- 
ity function of the globular clusters at 1'.5 < r < 
8' derived from the KPNO photometry, the lumi- 
nosity function of the background objects with the 
same color range as that of the globular clusters at 
9' < r < 10', and the net luminosity function for 
1'.5 < r < 8' that was derived after subtracting 
the background luminosity function from the lumi- 
nosity function of the globular clusters. The lumi- 
nosity functions were corrected for incompleteness 
using the completeness values derived in Section 2. 
We also plot the luminosity function of the glob- 
ular clusters at r < 1'.5 derived from the HST 
photometry. The HST luminosity function of the 
globular clusters at r < 1'.5 was derived by count- 
ing the globular clusters in a half circle of radius 
r < 1'.5 and doubhng the resulting counts. The 
KPNO luminosity function is similar to that of the 
HST luminosity function for the bright magnitude 
range Ti < 23.5, showing a peak at Ti w 23. 
For Ti > 23.5, the KPNO luminosity function 
falls below the HST luminosity function, indicat- 
ing that the KPNO photometry get s more incom- 



plete in the faint magnitude range. iLarsen et al 



(j200l[ ) derived the peak magnitude of the lumi- 
nosity function of the globular clusters in the cen- 
tral region of M60 from the HST/WFPC2 data: 
y(peak) = 23.58 ± 0.08, which corresponds to 
Ti = 23.38 adopting the color (C - Ti) = 1.7 
in the equation given in the previous section, and 
the foreground reddening Ay = 0.088. This value 
is consistent with the peak magnitude seen in the 
KPNO luminosity function. 

We derived an estimate of the total number of 



globular clusters by (a) counting the brighter half 
of the luminosity function which is almost com- 
plete and (b) doubling it after subtracting back- 
ground contribution, assuming that the luminos- 
ity function of the globular clusters is symmet- 
ric around the peak magnitude. The number of 
globular clusters brighter than the peak magni- 
tude (19 < Ti < 23.383) at r < 8' is 2048: 366 
at r < 1'.5 from the HST images, and 1682 at 
1'.5 < r < 8' from the KPNO images. We then 
derived the number of point sources with the same 
color as the globular clusters in the background 
region at 9' < r < 10' - 40. If this number is 
scaled according to the area ratio of the regions at 
r < 8' and at 9' < r < 10', it becomes 475. Part 
of these objects may be globular clusters, but we 
do not know what fraction of them. So we ar- 
bitrarily assume a half of these as globular clus- 
ters and assign the same value as an error. Then 
the total number of the globular clusters is esti- 
mated to be 2(2048 - 237.5) = 3621. We take 
N = 3600 ± 500 as the total number of globu- 
lar clusters in M60. This estimate is very sim i- 
lar to the value derived by iForbes et all (|2004h . 
3700 ± 900. From this value and the absolute 
magnitude of M60 {My = -22.44) we derive the 
specific frequency to be Sn = 3.8 ± 0.4, which i s 
similar to the value given bv lForbes et~al~l (|2004l ). 
Sn ~ 4.1 ± 1.0 , and substantially smal ler than the 
value given bv lAshman fc Zeof I (|l998h . Sn 6.7, 
based on inferior data. Our value, although on the 
small side for a cluster gE, is not atypical. 

4. Discussion 

4.1. Comparison vifith Previous Studies 

There is only one previous photometric study 
of the globular clusters of M60 cover ing a rela- 

12004,1 who 



tively large field, that of iForbes et a 



covere d about 90 square arcmin. IForbes et al. 
(l2004h showed that the color distribution of 995 
bright globular clusters with 20 < i < 23.6 and 
0.5 < {g — i) < 1.5 in M60 is bimodal, with peaks 
at {g-i) = 0.865±0.005 and 1.167±0.004. These 
peak colors correspond to (C — Ti) — 1.416 and 
1.989, respectively, using the transformation equa- 
tion between (g — i) and (C — Ti) given in Section 
2.5. These values are similar to those derived in 
this study, (C - Ti) = 1.37 and 1.87. Note that 
the color range between these two populations is 
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almost twice as large in (C — Ti) as {g — i), illus- 
trating the former 's superior metallicity sensitivity 
(jGeisler. Lee, fc Kim 1 996) . 

Forbes et al. ] (|2004 also presented the surface 



number density profiles of the globular clusters 
for the range of galactocentric distance 1'.3 — 3'. 8, 
showing that the RGCs with 0.5 < {g ~ i) < 0.9 
have a similar surface density distribution to that 
of stellar light, which is steeper than that of the 
BGCs with 1.1 < (g-i) < 1.5. They found that 
the surface number density profiles of the BGCs, 
RGCs and all globular clusters are fit well by 
power laws with slopes of -1.04±0.09, -1.73±0.06 
and —1.3 ± 0.05. Their slopes for the BGCs and 
all globular clusters are similar to our estimates, 
-1.128 ± 0.040 and -1.285 ± 0.034, but their 
slope of the RGCs is somewhat steeper than ours, 
-1.452 ±0.050. 



Forbes et al. 



(|2004l ) presented also the radial 
variation of the mean colors of the globular clus- 
ters for the range of galactocentric distance 0'.4 — 
6'. 3, showing that both the BGCs and RGCs show 
little radial color gradient, while the color of all 
globular clusters gets bluer with increasing radius. 
Our results on the radial gradient of colors shown 
in Figure 15 are consistent with these results. 

4.2. Metallicity Distribution 

Colors of old globular clusters with similar ages 
are determined primarily by the metallicity of the 
globular clusters. There are several studies for de- 
riving transformation relations between colors and 
[Fe/H] for globular clusters, which were based on 
the integrated photometry of the Galactic globular 
clusters, the combined data of globular clusters in 
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However, there are only a few studies on the 
transformation relation between (C — ri)o and 
[Fe/H] for glob ular clusters in the literature. 
Geisler fc Forte I f|1990. ) presented a linear rela- 



tion derived from the data for the Galactic glob- 
ular clusters: f Fe/H] = 2.35(C - T)o - 4.39; 
Harris &z Harrisi ( 2002[ ) derived a quadratic rela- 
tion based on the data for the Galactic globular 
clusters: [Fe/H]= -6. 037[1 - 0.82(C - rQo + 



(|2003l ) derived another quadratic relation based 
on the combined data including the globular clus- 
ters in M49 and M87 as well as the Galactic 
globular clusters: [Fe/H]= -0.75(C - ri)^ + 
4.438(C — ri)o — 5.64. However, the sequence 
for M87 globular clusters is not consistent with 
those of M49 globular clusters and Galactic glob- 
ular cl usters in the [Fe/H] versus (C — Ti)o dia- 
gram (jCohen. Blakeslee fc Cote 1 120031 ). while the 
sequences for the latter two are consistent. The 
[Fe/H] dat a for 38 glob u lar cl usters in M60 pre- 
sented by IPierce et al. ] (Hooi show too large a 
scatter to be used for this purpose. 

We derived new relations using the combined 
data for the globular clusters i n our Galaxy and 
M49. The [Fe/H] data are fro m [Harris I (Il996l) for 
Galactic globular clusters, and lCohen. Blakeslee fc Cote 
(|2003[ ) for M49 globular clusters. T h e (C - Ti) 



data are from Harris fc Canterna (Il977l) for 



0.162(C-ri)o ]; and ICohen. Blakeslee fc Cote 



Galac ti c globular clusters, and fromlGeisler. Lee, fc Kim 
(Il996l) : iLee. Kim fc Geisleil (Il998l) for M49 glob- 
ular clusters. We updated the (C — Ti) data for 
the Galactic globular clusters by adding our un- 
published data for two globular clusters, NGC 
6624 and NGC 6316, (C - Ti) = 2.75 ± 0.02, and 
2.096 ±0.02, respectively. We applied the redden- 
i ng co rrection using the values given in ,Harris1 
(|1996[ ). 

Figure 18 displays a [Fe/H] versus (C — Ti)o di- 
agram for these data. The sequence for the M49 
globular clusters (open squares) matches well that 
for the Galactic globular clusters, and covers as 
well the high [Fe/H] range where there are only 
a few data for Galactic globular clusters. The re- 
lation between [Fe/H] and (C — Ti)o appears ap- 
proximately linear over most of the range, with 
some possible curvature at both ends, especially 
at the metal-poor end. We tried to fit the com- 
bined data of our Galaxy and M49, after remov- 
ing outliers, using various equations, and found 
that a 3rd order polynomial and double linear 
relations give the best fits: [Fe/H]= 1.387(C - 
Ti)l - 6.698(C - Ti)l + 12.609(C - Ti)o - 9.379 
with reduced = 2-529, and [Fe/H] = (2.359 ± 
0.051)((C-Ti)o-1.46)-0.92for (C-Ti)o < 1.46, 
= (1.951 ± 0.044)((C - Ti)o - 1.46) - 0.92 for 
(C - ri)o > 1.46 with reduced = 5-090. These 
relations are similar in general to those given in the 
literature, showing some difference only around 
the low and high metallicity ends. But note that 
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the ICohen. Blakeslee fc Cote I (|2003l ) curve is sig- 
nificantly displaced from the other relations at vir- 
tually all metallicities. 

We derived [Fe/H] for M60 globular clus- 
ters from (C — Ti)o col ors u sing our relations 
and previous relation s (iGeisler fc Forte , ,1990t 



Harris fc l2002t ICohen. Blakeslee fc Cd~ 



20031 ). and displayed them in Figure 1 9. We 



plotte d also the [Fe/H] distribution from [Harris 



(Il996l) :or the Galactic globular clusters for com- 
parison in the same figure. Figure 19 shows 
several notable features as follows. First, the 
[Fe/H] distributions are very broad, covering from 
[Fe/H]R:i —2.4 to much higher than the solar 
value, reaching [Fe/H]w 0.8 and much higher 
than in the Galaxy. Second, all the [Fc/H] dis- 
tributions derived from various transformation re- 
lations show a dominant peak at [Fe/H]« —1.2 
dex. Third, none of the [Fc/H] distributions de- 
rived from various transformation relations look 
symmetric. Most of them show clearly a weaker 
component at [Fe/H]« —0.2 dex in addition to 
the domina nt peak at [Fe/H ]Ri 3 1-2 dex, al- 
though the [Harris fc HarrisI (|2002h calibration 
shows a dominant metal-rich peak. Thus they 
are all bimodal in metallicity. A typical photo- 
metric error of(T(C — Ti) = 0.1 leads to an error 
cr[Fe/H]= 0.20 ^ 0.24 when using the two trans- 
formation relations derived in this study. 

KMM tests of the data based on the 3rd order 
polynomial transformation show that the prob- 
ability that the [Fe/H] distribution is bimodal 
over unimodal is higher than 99.9 %. We deter- 
mined the parameters for best-fit double Gaussian 
curves to the metallicity distribution data for 1236 
globular clusters, using a maximum-likelihood 
method through the KMM mixture modeling 
routine (jAshman. Bird, fc: ZeDfl ll994V a metal- 
poor component with center at [Fe/H]= —1.18, 
a — 0.42, and N=779, and a metal-rich compo- 
nent with center at [Fe/H]= —0.27, a = 0.59, and 
N=467 for the 3rd order polynomial transforma- 
tion. Similar results are obtained for the double 
linear transformation: a metal-poor component 
with center at [Fe/H]= -1.22, a = 0.39, and 
N=747, and a metal-rich component with center 
at [Fe/H]= -0.18, a = 0.44, and N=489. Fourth, 
the metallicities of both peaks in M60 are 0.3 - 
0.4 dex larger, respectively, than those for Galactic 
globular clusters at [Fe/H]= —1.5 and —0.6. 



We compared [Fe/H] derived using the double 
linear r elation in this study with the estimates 
given bv lPierce et al. I ( 2006( i for 38 common glob- 
ular c lusters, as shown in Figure 20. iPierce et alTI 
(|2006l ) derived [Fe/H] of 38 globular clusters us- 
ing two methods from spectra: one using the 
simple stellar population (SSP) models and the 
other using the Brodie fc Huchra (BH) method 
(|Brodie fc Huchra Ill990( ). Fi gure 20 displays SSP 
and BH metallicities given bv lPierce et al. T (|2006l) 
versus [Fe/H] derived in this study. It is seen that 
two show g ood correlation but with some scatter. 
In addition, IPierce et al. 1 (l2006h 's [Fe/H] IS on av- 
erage 0.3 to 0.4 dex lower than our [Fe/H]. Linear 
fits to the data yield: SSP [Fe/H] (Pierce et al.) 
= 1.1 20 [Fc/H] (this study) -0.232 with rms=0.408, 
and BH [Fe/H] (Pierce et al.) = 0.805[Fe/H](this 
study) -0.558 with rms=0.315. The cause for this 
difference is not known. 

4.3. Comparison v^^ith M49 (NGC 4472) 

There are several gEs in Virgo and they are 
an excellent pool for studying the properties of 
the globular clusters in gE, because they belong 
to the same galaxy cluster, being located at the 
similar distance from us. M49 (NGC 4472) is the 
brightest gE in Virgo, located about 4 deg from 
the Virgo center. It is a prototypical example of 
a galaxy showing the bimodal color distribution 
of the globular clusters in gEs. Detailed anal- 
ysis of the Washingtong CTi photometry of the 
globular cluster s ystem in M49 based on si milar 
data was given bvlGeisler Lee, fc KimI (Il996f l and 
Lee. Kim fc Geisleij (|l998l ). Here we compared in 
detail the properties of the globular cluster sys- 
tems in M60 and M49, using very similar data and 
techniques. We have also limited ourselves to al- 
most identical globular cluster sample definitions: 
r = 1 - 7' and Ti = 19 - 23. The color boundary 
for M49 ((C - Ti) = 0.9 - 2.3) is 0.1 bluer than 
that for M60. 

Figure 21 displays the color distribution of the 
bright globular clusters in M60 and M49. It 
shows that the color distributions of the bright 
globular clusters in both galaxies are similarly bi- 
modal. We determined the parameters for best-fit 
double Gaussian curves to the color distribution 
data for M49, using a maximum-likelihood method 
through the KMM mixture modeling routine: a 
primary component with center at (C— Ti) = 1.30, 
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(7 = 0.13, and N=792, and a secondary compo- 
nent with center at (C-Ti) = 1.79, a = 0.21, and 
N=818. One difference between M60 and M49 is 
that the peak colors for M60 ((C - Ti) = 1.37 
and 1.87) are sUghtly redder than those for M49 
((C — Ti) = 1.30 and 1.79), whi ch is c onsistent 
with the finding bv lCouture et al.l ( 199lf ) based on 
BV photometry. The difference in the foreground 
reddening between M60 {E{B - V) ^ 0.026, 
E{C - Ti) = 0.051) and M49 {E( B -V)= 0.022 , 
E{C - Ti ) = 0.043) is negligible (|Schlegel et al 



Interestingly, previous studies have found signif- 



icant blue tilts in M87 and M60, 



Il998l ) . Therefore the colors of the globular clus- 
ters in M60 are considered to be on average intrin- 
sically about 0.1 mag redder than those in M49. 
If the globular clusters are of similar old ages, this 
implies that the M60 globular clusters are on aver- 
age about 0.2 dex more metal-rich than their M49 
counterparts. 

Recently IStrader et al. ( 2007^ found from the 
re-analysis of the spectr oscopic data for 47 bright 



globu lar clusters in M49 (|Cohen. Blakeslee fc Cote 



20031) that the metallicity distribution of these 
globular clusters is bimodal with two peaks at 
[m/H]~ —1.1 and 0.0. The peak metallicity for 
the metal-poor globular clusters in M49 is ^ 0.4 
dex higher than that for the metal-poor globular 
clusters in our Galaxy, and the peak metallicity 
for the metal-rich globular clusters in M49 is ^ 0.5 
dex higher than that for the metal-rich globular 
clusters in our Galaxy. These results are similar 
to those derived from the Washington photometry 
of the globular clusters in M49 and M60 in this 
study. 

The spatial distribution of the globular clus- 
ters in both galaxies share common features: (a) 
the RGCs are more centrally concentrated than 
the BGCs, (b) the globular cluster system is more 
extended than the stellar halo, and (c) the elon- 
gation of the RGC system is consistent with that 
of the stellar halo, while that of the BGC system 
is approximately circular, (d) the color gradient 
in the overall globular cluster systems are similar, 
as are the relative lack of color gradients in the 
individual RGC and BGC populations. 

Recent studies based on the analysis of HST 
ACS/WFC data for the globular clusters in ellip- 
tical galaxies found that the brighter the bright 
BGCs are, the redder their colors g et, which has 
been referred to as t he 'b l ue tilt' (IHarris et al. 
20061 : IStrader et al"! 120061 : iMieske et al. I l2006h . 



in M49 (S trader et al. 2006; Mieske et al 



but no ev idenc e 
20f)i ). 



All these results were based on the HST obser- 
vation of small fields most ly covering the cen- 



tral regions of the galaxies. iMieske et al. I (|2006l) 



found also that the blue tilts are seen in both in- 
ner regions (at r < 65") and "outer" regions (at 
r > 65") of bright galaxies, and that the slope for 
the globular clusters in the inner region is two-to- 
three times steeper in d{g — z)/dMz than that in 
the outer region. 

Figure 22 displays the color-magnitude dia- 
grams of M60 and M49 derived from the KPNO 
images for the outer regions a^J/^5_^£_^_7^^_We 
used the da ta for M49 from Geisler. Lee. &: Kim 
(Il996[ ) and iLee. Kim fc Geislerl (ll998V We de- 
termined the parameters for the best-fit double 
Gaussian curves describing the color distribution 
for given magnitude range, using a maximum- 
likelihood method through the KMM mixture 
modeling routine, and plotted the center values 
of the Gaussian components in Figure 22. 

In Figure 22 the blue tilt is clearly seen for the 
bright BGCs in M60, while it is not clearly seen 
for the bright RGCs. Only the bright BGCs with 
Ti < 22 in M60 show this blue tilt, which is con- 
sistent with the finding f o r the brightest cluster 
galaxies by iHarris et al. I (|2006l ) that the bright 
BGCs with Mi < -9.5 show the blue tilt, while 
the faint BGCs Mj > -9.5 do not. The blue tilt 
is also seen for the bright BGCs with Ti < 22 mag 
in M49, but it is much weaker than that for M60. 
The color change for the BGCs with 20 < Ti < 22 
is about 0.2 mag for M60, while it is about 0.1 
mag for M49. However, the brightest bin for M60 
is not separated clearly into two groups so that 
the difference in color change between M60 and 
M49 may not be significant. This is in contrast 
with the case of the inner region of M49 where 
no blue tilt was found in the HST ACS images 
(jStrader et al. II2006I : 



Mieske et ani2006f) . Inter- 



estingly the RGCs in M49 also show a similar tilt 
like the BGCs, which is different from the typical 
CMD. Therefore it is needed to investigate this 
problem with better data. 

Several mechan isms were discussed for the ori- 
gin of th e blue tilt (IHarris et al. l2006t Strader et al. 
.2006; M ieske et al l2006tlBekki. Yahagi. fc Forbes I 
[20071): accretion of globular clusters from low- 
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mass galaxies, the capture of field stars by individ- 
ual globular clusters, self-enrichment in globular 
clusters, contamination effect (due to super-star 
clusters, stripped nuclei, and ultra compact dwarf 
galaxies (UCDs)) , and stochas t ic effe cts. After 
examining these, iMieske et al. ( 2006 ) concluded 
that self-enrichment and field star capture are the 
most promising to explain the observational re- 
sults for the blue tih derived from the HST ACS 
data for a sample of 79 early-type galaxies. 

We discuss these two mechanisms in relation 
with the observational results for M60 and M49. 
If the blue tilt is mainly due to field star capture, 
it is expected that there is a radial variation of the 
blue tilt in the sense that it is steeper in the inner 
region. This is because there is a negative radial 
color gradient for halo stars at r < 2'. 5 as seen in 
the surface color profiles, and because more mas- 
sive globular clusters can collect more field stars 
that are redder in the inner region. However it is 
not possible to investigate this trend reliably with 
the current data for M60 and M49. 

On the other hand, the degree of self-enrichment 
is affected mainly by the mass of the stellar sys- 
tems. So the self-enrichment effect of the globular 
clusters is expected to depend on the mass of the 
globular clusters, but little on the galactocentric 
distance. This is consistent with the observational 
results for the outer regions of M60 and M49. 
However, the existence of self-enrichment is seen in 
few, if any, typical globular cl usters in our Galaxy 
and is a controversial subject ( Thoul et al. I (|2002l ) 
and references therein). It is also difficult to ex- 
plain the absence of the blue tilt in the inner region 
of M49 with self-enrichment. Further studies of 
globular clusters covering the full spatial extent of 
the globular cluster systems of gEs are needed to 
understand this problem. 

4.4. Implications for the Formation of 
Globular Clusters 

Key observational facts for globular clusters in 
M60 to be explained by any formation scenario 
can be summarized as follows. 1) The color distri- 
bution of the globular clusters is bimodal. 2) This 
indicates that the range of metallicity of the glob- 
ular clusters is very large from [Fe/H]< —2.0 to 
~ 0.5, and that the metallicity distribution is also 
bimodal, consisting of a strong metal poor compo- 
nent with peak [Fe/H]=— 1.2 and a weaker metal 



rich component with peak [Fe/H]Ri —0.2. 3) The 
spatial distribution - density, ellipticity and posi- 
tion angle - of the RGCs is similar to that of halo 
stars, while the spatial distribution of the BGCs 
is much more extended and circular than that of 
the stellar halo. 4) The bright BGCs show a blue 
tilt, while the bright RGCs do not. 

Several scenarios to explain the formation of 
globular clusters i n gEs have been propos ed: pri- 
mordial formation ( Peebles fc Dicke 1 ll968l) includ- 
i ng th e biased globular cluster formation bv'West] 
(|l993l ). gaseous merger (Ashm an fc Zepf , .1992r). 



in situ multiphase collapse mo del ([Forbes et al 



19971 : iKissler-Patig et al 



'1998'), and tidal strip- 
ping/accretion (jCote et al. 1998). Each of the 
scenarios has pros and cons, but none of them by 
itself can explain all the observational results. The 
current paradigm in the formation scenarios of 
the globul ar clusters a ddresses the following three 
processes '2003; 'Kravtsov fc Gnedin 

Brodie fc Strader 2006: Harris et al 



g tnree 
T I2OO5I: 



1. Formation of the first generation of globu- 
lar clusters: The metal-poor globular clus- 
ters were formed first a mong the stellar 



syste ms in the universe ([Peebles fc Dicke 



1968]). The mass of their progenitors must 



be much larger than the mass of the present 
globular clusters, possibly similar to that of 
dwarf galaxies (10^ — 10^ Mq), to explain 
the blue tilt via self enrichment. They were 
formed not only around massive galaxies, 
but also in the field throughout the universe 
where the gas density was sufficiently high. 
The latter could now be intracluster glob- 
ular clusters. However, the fact that the 
metal-poor globular clusters are found in 
and around the galaxies implies that they 
were formed preferentially around the high 
density peaks rather than in any r andom 



field i n the universe, as proposed by iWest 
([19931) . 



2. Formation of the second generation of globu- 
lar clusters: The metal rich globular clusters 
were formed together with the halo stars, 
both of which were formed in the interstel- 
lar medium that was enriched rapidly after 
the formation of the first generation of glob- 
ular clusters. This enrichment time scale 



should be smaller than about two Gyr ([Lee 
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20031 ). and the enrichment range should be 
large, from [Fe/H]< —1.2 to « —0.2, a range 
similar to that between the metal-poor and 
metal-rich Galactic globular cluster popula- 
tions but at a significantly higher metallic- 
ity. The enrichment process can be double- 
episodal or continuous, depending on the ex- 
istence of two peaks or a singl e broad peak in 



the m etallicity distribution (jBeaslev et al 



20021 ). The broad metallicity distribution 



and the high metallicity peak being much 
weaker than the low metallicity peak, if any- 
thing, indicate that the enrichment process 
might have been mainly continuous rather 
than double-episodal. However, the exis- 
tence of a minor high metallicity component 
indicates that there is a separate component 
of globular clusters that were formed in an- 
other episode. The color gradient of these 
globular clusters is very similar to that of the 
halo stars. These globular clusters in giant 
elliptical galaxies can be formed via gaseous 
merger of galaxies. If they are formed via 
gaseous merger, they must have been formed 
during the early phase of merging (within 
about two Gyr from the birth of their host 
galaxies). Later merging can make the host 
galaxi es grow mostly in mas s, not in metal- 
licity (|De Lucia et al. It is expected 
that the spatial distribution of these metal- 
rich globular clusters and halo stars is more 
centrally concentrated compared with that 
of the metal-poor globular clusters. The spa- 
tial distribution of the metal-rich globular 
clusters can be similar to or slightly more 
extended than that of the halo stars, depend- 
ing on the formation time scale and dissipa- 
tion effect_o|_sta£sand_gl^ clusters (see 
also lTamura et al. I (|2006bh ). 



3. Tidal capture/ Accretion of the existing glob- 
ular clusters: Massive galaxies can collect 
globular clusters in the nearby surroundings 
that may include both lower mass galaxies 
or the field. Globular clusters accreted af- 
ter about one or two Gyr from the first star 
formation can have any range of metallic- 
ity. However, most of the accreted globu- 
lar clusters are probably metal-poor rather 
than metal-rich, since low mass galaxies have 
generally more metal-poor globular clusters 



(jMiller 1120061 ). and they are more likely to be 
seen in the outer regions of massive galaxies. 



Observational results for the globular clusters 
in M60 show features of all three processes. The 
metal-poor globular clusters are involved with the 
first and third processes, while the metal-rich glob- 
ular clusters are related with the second process. 
Therefore, all three mechanisms may have played 
a role in the formation of the M60 globular clus- 
ters. 

5. Summary and Conclusion 

We have presented Washington CTi photome- 
try of the globular clusters in M60 as well as that 
of the star clusters in NGC 4647, a companion 
spiral galaxy, covering a 16'. 4 x 16'. 4 field. We 
have analyzed various photometric properties of 
the globular clusters in M60. Primary results are 
summarized as follows. 

1. The color-magnitude diagram reveals a sig- 
nificant population of globular clusters in 
M60, and a large number of young luminous 
clusters in NGC 4647. 

2. The color distribution of the globular clus- 
ters in M60 is clearly bimodal, with a blue 
peak at (C — Ti) = 1.37, and a red peak at 
at (C-Ti) = 1.87. 

3. We derived two new transformation relations 
between the (C — Ti)o color and [Fe/H] us- 
ing the data for the globular clusters in our 
Galaxy and M49: a 3rd order polynomial re- 
lation, [Fe/H]= 1.387(C - Ti)l - 6.698(C - 
Ti)2 + 12.609(C - Ti)o - 9.379, and a dou- 
ble linear relation, [Fe/H] = 2.359((C - 
Ti)o - 1.46) - 0.92 for (C - Ti)o < 1.46, 
and = 1.951((C - Ti)o - 1.46) - 0.92 for 
(C — ri)o > 1-46. Using these relations 
we derived the metallicity distribution of the 
globular clusters in M60, which is bimodal: a 
dominant metal-poor component with center 
at [Fe/H]w —1-2, and a weaker metal-rich 
component with center at [Fe/H]« —0.2. 

4. The radial number density profile of the 
globular clusters is more extended than that 
of the stellar halo. The radial profiles for 
the outer region at 1' < r < 7' are fit 
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approximately equally well by the deVau- 
couleurs law and a power law. The radial 
number density profile of the BGC is more 
extended that that of the RGCs. The ra- 
dial profiles for the central region show flat- 
tening, and the radial profiles for r < 2' 
is well fit by the King model. The core 
radii derived for the fix concentration value 
of c = 2.5 are Ire = 0.85' for all globu- 
lar clusters, rc = 1.06' for the BGCs, and 
Tc = 0.72' for the RGCs. Thus the core ra- 
dius for the RGCs is much smaller than that 
for the BGCs. 

5. The surface number density maps of the 
globular clusters show that the spatial dis- 
tribution of the BGCs is roughly circular, 
while that of the RGCs is elongated simi- 
larly to that of the stellar halo. 

6. The mean color of the bright BGCs gets red- 
der as they get brighter both in the inner and 
outer region of M60. This blue tilt is seen 
also in the outer region of M49, the brightest 
Virgo galaxy. 

7. We estimated the total number of the glob- 
ular clusters in M60 to be 3600 ± 500, and 
the specific frequency to be Sn = 3.8 ± 0.4. 
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Table 1 

Basic information of M60 (NGC 4649) 



Parameter 


Values 


References 


RA(2000), Dec(2000) 


121i43m39.66s, +11° 33' 9".4 


1 


Effective radius, Reff 


97 arcsec (C), 110 arcsec (Ti) 


2 


Effective ellipticity, eg// 


0.21 (C,Ti ) 


2 


P.A. 


106 deg (C), 105 deg (Ti) 


2 


Standard radius, R25 


242 arcsec (C) 


2 


Standard ellipticity, 625 


0.224 (C) 


2 


P.A. (i?25) 


108 deg (C) 


2 


Total magnitudes 


= 8.84 ± 0.05, = 9.81 ± 0.05 


3 


X-ray luminosity 


LogiLx) = 41.16 [erg s~i] 


4 


Systemic radial velocity, Vp 


1117 ±6 km s-i 


1 


Foreground reddening 


E{B -V)^ 0.026 


5 


Distance 


d=17.30 Mpc ((m - M)o = 31.19 ± 0.07) 


6 



References. 



(1 ) NA SA Ext ragalactic Database; (2) 



de Vaucouleurs et al. ] (119911): (4 ) lO'Sullivan. Forbes, fc PonmanI 
Schleeel et al. , C1998i); (6) iMei et al.1 (|2007l) . 



This study; 
( 200lh : 



(3) 
(5) 
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Table 2 
Observation log 



Target 


Filter 


T(exp) 


Airmass 


Seeing 


Date(UT) 


M60 


C 


100 


1.4 


1".43 


1997 Apr 9 


M60 


C 


1500 


1.5 


1".47 


1997 Apr 9 


M60 


R 


60 


1.3 


1".43 


1997 Apr 9 


M60 


R 


1000 


1.4 


1".35 


1997 Apr 9 


M60 


C 


3 X 1500 


1.2 


1".09 - 1".24 


1997 Apr 10 


M60 


R 


2 X 1000 


1.4 


1.34-1".44 


1997 Apr 10 



Table 3 

cti photometry of the measured point sources with a{c — ti) < 0.3 in mso" 



ID 


X[px]^ 


Y[px]'' 


RA[dcg] 


Dcc[dcg] 


Ti 




[C - Ti) 


<j{C - Ti) 


2 


182.38 


457.14 


190.800629 


11.479574 


17.714 


0.039 


3.265 


0.039 


4 


991.59 


1164.33 


190.908737 


11.571802 


17.898 


0.033 


2.386 


0.038 


5 


665.44 


1204.13 


190.865173 


11.577062 


17.925 


0.025 


1.095 


0.031 


6 


144.16 


112.40 


190.795517 


11.434722 


18.048 


0.063 


3.753 


0.066 


7 


862.64 


923.39 


190.891418 


11.540276 


18.060 


0.014 


2.204 


0.028 



'^The complete version of this table is in the electronic edition of the Journal. The printed 
edition contains only a sample. 

'^X and Y increases toward east and north, respectively. A pixel corresponds to 0.47 arcsec. 
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Table 4 

Number density profiles of globular clusters with 19 < Ti < 23 in M60 



r[arcmin] 


cr(All GC)'^ 


ct(BGC)" 


(t(RGC)'^ 


0, 


.083 


115.870 ±37.203 


45.474 ±23.529 


70.394 ±28.818 


0, 


.292 


85.695 ± 15.463 


33.950 ±9.856 


51.742 ± 11.915 


0, 


.542 


78.312 ± 12.663 


34.455 ±8.495 


43.855 ±9.391 


0, 


.833 


47.184 ±6.602 


20.670 ± 4.451 


26.512 ±4.876 


1, 


.167 


29.771 ± 4.756 


14.189 ±3.363 


15.580 ±3.363 


1, 


.750 


20.160 ±2.082 


10.668 ± 1.556 


9.490 ± 1.383 


2, 


.250 


14.561 ± 1.626 


6.988 ± 1.178 


7.571 ± 1.121 


2, 


.750 


10.288 ± 1.258 


5.496 ± 0.964 


4.789 ± 0.808 


3, 


.250 


9.890 ± 1.112 


5.945 ±0.894 


3.942 ±0.662 


3, 


.750 


5.714 ±0.799 


2.586 ±0.596 


3.126 ±0.532 


4, 


.250 


5.339 ± 0.730 


3.434 ±0.613 


1.903 ±0.396 


4, 


.750 


5.059 ±0.677 


3.241 ±0.569 


1.817 ±0.367 


5, 


.250 


3.921 ±0.588 


2.780 ±0.514 


1.140 ±0.284 


5, 


.750 


3.983 ± 0.565 


2.846 ± 0.495 


1.135 ±0.271 


6, 


.250 


1.990 ±0.438 


1.367 ±0.388 


0.620 ±0.204 


6, 


.750 


2.940 ± 0.472 


1.906 ±0.406 


1.032 ±0.241 


7, 


.250 


2.086 ±0.412 


1.620 ±0.375 


0.465 ± 0.170 


7, 


.750 


1.174 ±0.355 


0.425 ±0.293 


0.747 ±0.201 


8, 


.250 


1.040 ±0.411 


0.632 ±0.365 


0.405 ±0.190 


8, 


.750 


1.139 ±0.493 


0.831 ± 0.449 


0.306 ± 0.204 



"^Number per square arcmin after background subraction. 
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Fig. 1. — A grayscale map of the Washington 
Ti image of M60. The size of the field of view 
is 16'. 4 X 16'. 4. North is up and east to the left. 
The small spiral galaxy in the north-west (ARA= 
-109".9 and ADec= 107". 3) of M60 is NGC 4647. 
The HST/WFPC2 fields are also marked by boxes. 



Fig. 2. — Mean photometric errors of Ti and 
(C — Ti) vs. Ti magnitudes of the point sources 
measured in the long exposure images of M60. In 
the low panel, squares, circles, and crosses repre- 
sent the errors for all sources, blue globular clus- 
ters (1.0 < (C-Ti) < 1.7), and red globular clus- 
ters (1.7 < (C - Ti) < 2.4), respectively. 
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Fig. 3.— Completeness for KPNO C (a) and Ti 
(c): thick solid line for the entire region, short dash 
line for 1' < r < 2', long dash line for 2' < r < 3', 
dot-short dash line for 3' < r < 4', dot-long dash 
line for 4' < r < 5', and short dash- long dash line 
for 5' < r < 6'. Completeness for HST V (b) and 
/ (d): thick solid line for the entire region, short 
dash line for r < 10", long dash line for 10" < 
r < 25", dot-short dash line for 25" < r < 40", 
dot-long dash hne for 40" < r < 60", and short 
dash- long dash line for 60" < r < 80". 
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Fig. 4.— (a) iForbes et all (|2004f )'s {g - i) vs 
(C— Ti) of the measured obje cts with Ti < 2 2 maj 
in M60 in common between iForbes et al. I ( 200 
and this study. The solid lines represent multi- 
linear fits to the data, (b) HST {V-I) vs. (C-Ti) 
of the measured objects with V <22 mag in M60. 
The solid line represents a linear fit to the data, 
and the dashed line represents the fit derived from 
the photometry of M49 (NGC 4472) by Lee & Kim 
(2000). 



23 



20 



24 



S6 



_,, 1 11 11 1 11 11 1 11 11 1_ 

- (a) HST(r<1.5') - 


18 










20 










22 










24 


1 


2 


3 







Fig. 5. — Color-magnitude diagrams of the point 
sources in M60. The dashed line represents an 
approximate lower boundary of the KPNO pho- 
tometry. The boxes represent the regions for se- 
lecting BGCs and RGCs with 19 < Ti < 23 
mag. (a) r < 1'.5 in the HST WFPC2 images. 
(C — ri)(-v_/) is the color converted from — /). 

(b) 1'.5 < r < 7' except for the circular region 
with radius of one arcmin centered on NGC 4647. 

(c) 7' < r < 9'. (d) Background (r > 9')- 
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Fig. 6. — Ti-{C — Ti) color-magnitude diagrams 
of the point sources at r < 0'.5 (filled circles) and 
0'.5 < r < 1' (open circles) from the center of 
NGC 4647 (upper panel), and the point sources in 
the "control field" with the same area at the same 
galactocentric distance from M60 as NGC 4647, 
but in the opposite direction (ARA = 109.9" and 
ADec = —107.3") (lower panel). In the upper 
panel, filled squares and open squares represent 
slightly extended sources at r < 0'.5 and 0'.5 < 
r < 1.0', respectively. Most of these sources are 
considered to be star clusters associated with NGC 
4647. The boxes are the same as in Fig. 4. 
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Fig. 7. — (C — Ti) color distribution of the globu- 
lar clusters with 19 < Ti < 23 mag at 1' < r < 7' 
from M60. The thick solid line represents the dou- 
ble Gaussian fit to the data, each component of 
which is plotted by the dashed lines. The hatched 
histogram represents the color distribution of the 
background objects at 9' < r < 10'. 
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Fig. 8. — Radial variation of the (C — Ti) color 
distribution of the globular clusters with 19 < 
Ti < 23 mag in M60 after background subtrac- 
tion. (Left panel) KPNO photometry. (Right 
panel) HST/WFPC2 photometry. 
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Fig. 9. — Radial variation of the number ratio of 
the BGCs and the RGCs with 19 < Ti < 23 mag 
in M60. The solid line represents a linear fit to the 
data. 
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Fig. 10. — Surface photometry of M60. 
represents galactocentric distance along the ma- 
jor axis, (a) C and Ti surface brightness mag- 
nitudes vs. rl^a-j- The solid hnes represent fits 



with law for 4" < r < 410". (b) Surface 
(C — Ti) color vs. rmaj- The solid line along 
the data represents a fit to the data: (C — Ti) = 
-0.162 ± O.Oniogr + 2.160 for 3" < r < 151" 
(U — R ) and {B — R) colors given bv lPeletier et al. I 
are also plotted for comparison, (c) Ellip- 



ticity vs. rjnaj- (d) Position angle (PA) vs. 
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10 5 -5 -1010 5 -5 -10 

4RA [arcmin] A[!A [arcrnin] 

Fig. 11. — Spatial distribution of the globular 
clusters with 19 < Ti < 23 mag. The radii of the 
concentric circles are from 1 to 9 arcmin increas- 
ing by one arcmin, and a small circle of radius 
1 arcmin to the north-west of M60 is for NGC 
4647. (a) All globular clusters, (b) Foreground 
stars with (C - Ti) > 2.5 and Ti < 23 mag. (c) 
BGCs (open circles) and very blue objects with 
—0.5 < (C — Ti) < 0.5 (crosses, mostly in the 
region of NGC 4647). (d) RGCs. 
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Fig. 12. — Number density contours of the globu- 
lar clusters with 19 < Ti < 23 mag in comparison 
with the stellar halo of JVI60. The size of the field 
is 16'. 4 X 16'. 4. Contour levels are the same in 
(a), (c) and (d): 3.323, 6.647, 9.970, 13.294, and 
16.617 globular clusters per arcmin^. The cross 
represents the center of JVI60. (a) All globular clus- 
ters, (b) A gray scale map of the Ti image of M60, 
with isophotal contours. Circles represent bright 
foreground stars, (c) BGCs. (d) RGCs. 
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r [arcmin] 

Fig. 13. — Number density profiles of the globu- 
lar clusters with 19 < Ti < 23 mag (open circles 
for all globular clusters, open triangles for BGCs, 
and open squares for RGCs). Thick-dashed lines 
represent the C surface brightness of the stellar 
halo of M60, converted using (26 — ^c')/2-5. Solid 
lines represent fits to the three different datasets 
for 0'.5 < r < 7'. (a) Number density vs. r^/'^ 
[arcmin]. (b) Number density vs. r [arcmin]. 
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Fig. 14. — Number density profiles of the globular 
clusters with King model fits (open circles for all 
globular clusters, open triangles for BGCs, and 
open squares for RGCs). Curved lines represent 
King model fits to the three different datasets for 
r < 2' (marked by the vertical dotted line), (a) 
Error weighted fit. (b) Equal-weighted fit. (c) Fit 
with the concentration parameter c = 2.5. 
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r[arcmin] 

Fig. 15. — (C — Ti) color vs. galactocentric dis- 
tance for the globular clusters with 19 < Ti < 23 
mag (dots) in comparison with the (C— Ti) surface 
color profile of stellar halo (curved line) in M60. 
Larger filled and open symbols represent, respec- 
tively, the mean and median colors in a given ra- 
dial bin, and solid lines are linear fits to the mean 
colors. Triangles, circles and squares represent, re- 
spectively, RGCs, all globular cluster, and BGCs. 
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Fig. 16. — Ti (left panel) and C (right panel) 
magnitudes vs. galactocentric distance for the 
globular clusters with 19 < Ti < 23 mag in M60. 
(a) All globular clusters, (b) BGCs. (c) RGCs. 
Filled circles represent the mean magnitudes in a 
given radial bin, and solid lines are linear fits to 
the mean magnitudes. 
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Fig. 17. — Luminosity function of the globular 
clusters in M60: dashed line for the raw luminosity 
function of the globular clusters at 1'.5 < r < 8' 
derived from the KPNO data, open squares for 
the luminosity function corrected for incomplete- 
ness, histogram for the luminosity function derived 
from the HST data for r < 1'.5, and cross for 
the luminosity function for the background ob- 
jects with the same color as the globular clus- 
ters at 9' < r < 10' normalized to the area for 
1'.5 < r < 8'. Filled circles represent the net lu- 
minosity function for 1'.5 < r < 8' after subtract- 
ing the contribution due to the background ob- 
jects. The HST luminosity function corrected for 
the completeness was scaled arbitrarily to match 
the KPNO luminosity function. 
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Fig. 18.— [Fe/H] versus (C - Ti)o for the 
globular clusters in our Galaxy (filled circles 
for E(B - V^) < 0.4, and filled triangles for 
E{B - V) > 0.4) and M49 (open squares). 
Lines represent the transformat ion relations be- 
tween [Fe/H] and { C - Ti)^ by iGeisler fc Forte I 
(|l990t ) (dotted l ine). [Harris fc Harris! (|2002|) (short 
dashed curve), ICohen. Blakeslee fc Cote I ( 2003h 
(long dashed curve), and double linear relations 
(thick dashed-dotted line) and 3rd order poly- 
nomial (thick solid curve) derived in this study. 
Large open circles indicate the globular clusters 
that were not used for the derivation of the trans- 
formation relations in this study. Mean errorbars 
are plotted in the lower right region. 
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Fig. 19.— Metallicity [Fe/H] distribution of the 
globular clusters in M60 (19 < Ti < 23 and 1' < 
r < 7') (circles), in comparison with the Galactic 
globular clusters (histogram). [Fe/H] was derived 
from the (C — Ti)o color using the following trans- 
formation relations: 3rd order polynomial equa- 
tion (circles in (a)), doubl e linear equation ( circles 



2002 ) 



2003^ 



in (b)) in this s tudy, and iHarris &: Harri; 
(circles in (c)), ICohen. Blakeslee fc Cote . 
(triangles in (c)), and lCeisler fc Forte 1 11990) (cir- 
cles in (d)). In (a) and (b) smooth solid lines 
and dashed lines represent the Gaussian fits to the 
data. 
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[Fe/H](This study) 
Fig. 20. — Compariso n of sp ectroscopic [Fe/H] 



given by IPierce et al. ( 20061 ) and photometric 



[Fe/H] derived using the double linear relation in 
this study for 38 common globular clusters in M60. 
Filled circles and open squares represent, respec- 
tively, [Fe/H] derived using the simple stellar pop- 
ulation (SS P) models and Brodi e & Huchra (BH) 
method by IPierce et al. ( 20061 ). Solid line and 
dashed line represent linear fits to the data for 
SSP [Fe/H] and BH [Fc/H], respectively, and dot- 
ted line represents one-to-one correspondence. 
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Fig. 21. — Comparison of the color distribution of 
the globular clusters with 19 < Ti < 23 mag and 
1' < r < 7' in M60 (this study, upper panel) and 
M 49 (Geisler, Lee, & Kim 1996, lower panel). 
Solid lines represent double Gaussian fits, while 
the dashed lines represent the individual Gaussian 
components. 
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Fig. 22. — Comparison of the radial variation 
of the color magnitude diagrams of the globular 
clusters in M60 (this study - left panel) and M 
49 (Geisler, Lee, & Kim 1996 - right panel). The 
data are plotted for the objects at 1'.5 < r < 7' of 
both galaxies. Dotted boxes represent the bound- 
ary for the globular clusters, and vertical dashed 
lines represent the peak values for the BGCs and 
RGCs. Circles and squares represent the Gaussian 
centers for BGCs and RGCs, respectively, derived 
for the magnitude bins using the KMM mixture 
model. 
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